Hierarchical clustering and the baryon distribution in galaxy clusters 



Eric R. Tittley and H. M. P. Couchman 

University of Western Ontario 

Department of Physics and Astronomy London, Ontario, N6A 3K7, Canada 

The baryon fraction of galaxy clusters in numerical simulations is found to be dependant on 
the cluster formation method. In all cases, the gas is anti-biased compared with the dark matter. 
However, clusters formed hierarchically are found to be more depleted in baryons than clusters 
formed non-hierarchically. There is a depletion of 10 to 15% for hierarchically formed clusters 
while the depletion is less than 10% for those clusters formed non-hierarchically. This difference is 
dependent on the mass of the clusters. The mean baryon enrichment profile for the hierarchically 
formed clusters shows an appreciable baryon enhancement around the virial radius not seen in the 
clusters formed without substructure. If this phenomenon applies to real clusters, it implies that 
determinations of the baryon fractions in clusters of galaxies require data extending beyond the 
virial radius of the clusters in order to achieve an unbiased value. 
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1. INTRODUCTION 

Since the early hydrodynamic simulations of galaxy clusters (Evrard 1990), the global distri- 
bution of the baryons has been of interest. In particular, the local enhancement or deficit of the 
baryons compared with the dark matter has profound implications for inferences of the universal 
baryon fraction as parametrized by ^li). Recall that is the ratio of the mass-energy density of the 
universe to the critical density required to close the universe. The ratio, Qf,, is the contribution of 
the baryonic component to 17, while 0„j is the contribution of all mass. The fraction, Qf,, is con- 
strained by primordial nucleosynthesis calculations which are sensitive to the cosmological model. 
Currently, estimates vary between a reported low of Cli, = 0.013 it 0.003^"^ (White and Fabian 
1995) to Ub = (0.020 ib0.002)/i~^ (Bludman 1998). Since galaxy clusters are the largest objects in 
the universe for which one can observe the baryon content as well as derive the total mass, they 
provide the most unbiased samples from which to calculate the baryon fraction of the universe. 
The values for fife and the baryon fraction, ft, = ^'>/n^, derived from observations of clusters, may 
be combined to derive 0^, the total contribution of matter to O. Observations of this sort find 
baryon fractions of 10 to 22% (White and Fabian 1995; White, Jones, and Forman 1997). This 
implies < 0.9 for = 0.1 and h = 0.5 while suggesting it is probably closer to 0.3 if we take 
= 0.15 and h = 0.65. However, the derivation of from the observations of clusters assumes the 
dark matter and hot gas are distributed in constant proportions within the cluster. Consequently, 
an understanding of the concentration factor of the gas in clusters is required since it represents 
the degree of biasing. As well, this biasing can have consequences in regards to constraining the 
deceleration parameter, qo (Rines, Forman, Pen, and Jones 1998). 

Numerical simulations generally agree that the gas is anti-biased in clusters of galaxies (Evrard 
1990; Thomas and Couchman 1992; Cen and Ostriker 1993; Kang, Cen, Ostriker, and Ryu 1994; 
Metzler and Evrard 1994; Pearce, Thomas, and Couchman 1994; Navarro, Frenk, and White 
1995; Anninos and Norman 1996; Lubin, Cen, Bahcall, and Ostriker 1996; Pildis, Evrard, and 
Bergman 1996). Results have been produced to the contrary (Owen and Villumsen 1997) using 
two-dimensional codes. Pearce et al. (1994) explains the anti-bias as a result of the merging pro- 
cess in which gas is shocked, permanently removing energy from the dark matter component and 
passing it to the gas. 

However, there is an uncertainty in the actual amount of bias. The bias between the universal 
baryon fraction and that found in clusters cannot be measured directly but must be inferred from 
numerical simulations. The numerical simulations must assume a model for the universe. Since the 

actual cosmogony of the universe is unknown in the details, variation among baryon fractions found 
in different model universes adds uncertainty. A better understanding of how this bias depends on 
cosmological models would assist in understanding the actual baryon fraction of the universe. 

In hierarchical clustering, the largest structures forming at a given time do so via the amalga- 
mation of many smaller structures which have formed at an earlier time. This is owing to the form 
of the initial density perturbation spectrum in which small-scale perturbations have higher initial 
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amplitudes than large-scale. In contrast, non- hierarchical clustering involves structure formation 
from the collapse of large structures with smooth density distributions. The results of numerical 
simulations compared with observations support the theory that we live in a universe in which 
structure is formed hierarchically. The degree to which the hierarchical nature affects the baryon 
distribution in galaxy clusters is not entirely clear. In this paper, a comparison of the baryon 
fraction bias will be made using numerical simulations of galaxy clusters formed hierarchically and 
non-hicrarchically. By Tising these extreme cases, the significance of hierarchical clustering itself 
will be determined. The analysis will concentrate not on the properties of individual clusters, but 
on the global mean properties of scaled quantities. 

The layout of the paper is as follows. In Sec. 2, the simulations are described. A brief 
description of the analysis methods are given in Sec. 3. The effect of hierarchical clustering on the 
baryon fraction is examined in Sec. 4. The results are discussed in Sec. 5. 
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2. SIMULATIONS 

The simulations used for this study are the same as described in Tittley and Couchman (1999). 
They comprise a set of five high resolution simulations in which hierarchical clustering was reduced 
or eliminated in four of the simulations via smoothing of the density perturbations in the initial 
conditions file. 

The initial, unsmoothed density distribution was scale-invariant, with a spectral index of n = 
— 1, such that the density fluctuation power-spectrum had the form of the power law, P{k) oc A;". 
It was normalised to set the variance on the scale of 8^"^ Mpc to ag = 0.935. 

The smoothing was done using a top-hat method for two of the simulations and a low-bandpass 
filter for the other two. For the top-hat method, smoothing over radii of 3h~^ Mpc and 7h~^ Mpc 
was done. The frequency cutoffs were 2T:{7h~^ Mpc) and 27r(14/i~^ Mpc). Subsequent references 
to these simulations will omit the and 27r factors. The Hubble constant parameter, h = 0.65, 
was adopted for these simulations. The details of the simulations are given in Table 1. 

The simulations contained 2 x 64^ particles, half of which were collisionless particles represen- 
tative of the dark matter component while the other half were collisional which represented a gas 
phase. The particles had masses of 9.27 x 10^° M© for the dark matter and 1.03 x 10^° M© for 
the gas. The total mass density is sufEcient to produce a flat cosmology, without a cosmological 
constant. 

The matter was evolved in a box with periodic boundary conditions on a scale of 40^~^ Mpc in 
co-moving coordinates. The N-body AP'^M — SPH code, HYDRA (Couchman, Thomas, and Pearcc 
1995) was used for all simulations. Two-body interactions were reduced by the use of gravitational 
softening. This was provided by the softening parameter, e = 20/i^^ kpc. Cooling was neglected. 
The cooling time for the bulk of the cluster gas is estimated to be well over the age of the clusters. 
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3. ANALYSIS 

The same sample of clusters as selected in Tittley and Couchman (1999) was used here. The 
details of the sets of clusters for each simulation are given in Table 2. Clusters were selected by a 
tomographic deprojection method. 

For each cluster, masses and radii were calculated at overdensities of 200 and 500. The radius 
for an overdensity of 200, -R20O) is the radius of a cluster interior to which the mean density is 
200 times the background density, which for these simulations is the critical density, pc- A similar 
definition holds for i?5oo- The mass interior to these radii are denoted, M200 and M500, respectively. 
The overdensity of 200 was chosen since it is on the conservative side of 178, the overdensity at 
which a spherically symmetric sphere of matter virialises, according to analytic models of top-hat 
collapses. The masses of both the dark matter, M^^rk, and gas, Mgas, were found at this virial 
radius, along with the total mass, M2oo- 

A lower limit to the size of the clusters was set by the requirement that each cluster, within 
the overdensity radius of R500, have at least 300 gas particles and 300 dark matter particles. This 
ensures the densities are calculated correctly (Tittley and Couchman 1999). 

Profiles of the dark matter and gas densities were calculated for each of the clusters. This was 
done by summing the mass contributions of the particles falling in exponentially separated radial 
bins centred on the clusters. The profiles were scaled by the approximate virial radius, -R200- 

Mean density profiles for each simulations were then calculated from the profiles of the indi- 
vidual clusters in each simulation. For the calculation of the mean density profile, the individual 
cluster profiles were weighted by their respective mass, M2oo- This weighting removes the bias 
towards the more plentiful low-mass clusters. However, for the density profiles, this weighting was 
found to have only a small effect, demonstrating the lack of dependency on mass for the profiles. 
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4. BARYON FRACTION 

4.1. The concentration parameter, T 

If wc define T200 cis the gas fraction in a spherical region of ovcrdensity 200 normahsed to the 
cosmic value, then we have (White, Navarro, Evrard, and Prenk 1993) 

Mgas + Mdark 

Similarly, we may define T(r) as the normalised gas fraction in a shell of radius r. Here, r is 
replaced by the dimensionless parameter r/Rg since radial profiles of both the dark matter and gas 
densities are characterised by the virial radius. 



4.2. Variation of T on cluster mass 

The values of T200, when plotted versus -M200 (Fig- 1), support the values found in Evrard 
(1997) and Eke, Navarro, and Frenk (1998) of T200 = 0.85 to 0.90 if the sample is restricted to the 
largest clusters. 

For the lower mass clusters, the clusters in the unsmoothed model have values ranging down 
to 0.75. There is no trend with mass for this model other than an increase in dispersion. The mean 
for the hierarchically formed clusters is 0.84 with the standard deviation varying from ±0.04 for 
clusters less massive than 10"*^^ Mq to ±0.01 for clusters more massive than 10^^ Mq . 

For the smoothed models, the trend is for T200 to decrease with mass with little change in the 
dispersion. For clusters more massive than 10^^ Mq , T = 0.87 to 0.97. Below 10^^ Mq , this drops 
down to the range 0.82 to 0.92. 



4.3. Baryon concentration profiles 

Since it is found that both the dark matter and gas density profiles scale with the virial radius, 
it stands to reason that the radial profile of the baryon concentration should also scale with radius. 
Numerical simulations of galaxy clusters using a 1-D Lagrangian code indicate that this is not true. 
Knight and Ponman (1997) reports that fh varies by a factor of about two at a given scaled radius 
between clusters that differ in mass by two orders of magnitude, with the less massive cluster having 
the smaller baryon fraction. This trend is consistent with the variation in T200 with the mass of 
the cluster seen for the clusters formed from smoothed initial conditions. However, since this factor 
is comparable to the cluster to cluster variation, it is justified once again to look at the mean radial 
profile of all the clusters in a model. 
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The baryon fraction profiles vary considerably for those clusters formed in the simulation with 
hierarchical clustering from those formed otherwise (Fig. 2, top) particularly around the virial 
radius. The baryon fraction is enhanced in a much deeper region (0.7 < r/i?200 < 10) in the 
hierarchical case than in the case with smoothed initial conditions (1 < r/i?200 < 4) and the 
enhancement is almost 4 times greater. In both cases, the peak in the baryon enrichment occurs 
just beyond the virial radius. Towards the more central regions, the enrichment in the hierarchical 
clusters drops off slightly more steeply. 

The discrepancies are not so large for the cumulative profile (Fig. 2, bottom). The clusters in 
the hierarchical model are generally more depleted in baryons than those in the smoothed models, 
in agreement with the results in Sec. 4.2. The mean profile for the unsmoothed model compares well 
those of the simulated clusters discussed in Eke et al. (1998). Those are simulated for a low-density 

universe. 

It must be noted that there arc large deviations in these profiles among the many clusters 
in the hierarchical scenario. These mean profiles represent trends. For the clusters formed non- 
hierarchically, however, the baryon profiles are fairly consistent among clusters. 

Interior to 1.5i?200) the mean T profile for the unsmoothed model may be fit by T(r) = 
0.67ib0.03+(0.52±0.05) ViJsoo ) non-cumulative. For the mean of the profiles for the non-hierarchically 
formed clusters, the fit takes a shallower form, T(r) = 0.76 ± 0.02 -|- (0.30 ± 0.07)7fl2„„. 
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5. DISCUSSION 

Using a scries of simulations in which clusters were formed either hierarchically and non- 
hierarchically, the relevance of hierarchical structure formation to the baryon distribution in clusters 
was explored. Non-hierarchical clustering was achieved by smoothing the initial density distribu- 
tion. 

It is found that the baryons are anti-biased in clusters of galaxies and this biasing is dependent 
on the presence of smoothing in the initial conditions. The bias parameter T ~ 0.85 for the 
hierarchically formed clusters. For the clusters formed from smoothed initial conditions, T ~ 0.92, 
indicating less of an anti-bias. For the non-hierarchically formed clusters, the normalised baryon 
fraction is dependent on the mass of the cluster, as well, which is not observed for the hierarchically 
formed clusters. 

In all cases, the baryon fraction was found to increase from the centre of the clusters, outward. 
Interior to 0.1i?200) the normalised factor spans T = 0.5 to 0.7. We determined that the baryon 
fraction profile, as parametrized by the normalised fraction, is steeper for hierarchically formed 
clusters. For these, it was found that oc (0.52zb0.05)r interior to 1.5i?2oo- For non-hierarchically 
formed clusters, the profiles were shallower, regardless of the degree of initial smoothing, with a 
mean dependency of oc (0.30 it 0.07)r. The baryon fractions peak between 1 to 2 virial radii. For 
the clusters formed hierarchically, the baryons 'pile up' in a deep region spanning 0.7 < r/i?200 < 10. 
This implies that measurements of the baryon fraction of clusters are sensitive to the radius outward 
to which the baryon content is integrated. 

Baryon fraction profiles determined from observations of clusters via the deprojection of the 
gas density (White et al. 1997) support the result that the baryon fractions increase with increasing 
radius (Markevitch and Vikhlinin 1997). However, the profiles are often shallower than those found 
here. White and Fabian (1995) finds from a sample of 19 clusters that the cumulative measurement 
of fb increases with radius following roughly the fb = 0.06i2 with R in Mpc . Similarly, White et al. 
(1997) finds a fit of fb = (0.12 ± 0.04)i2°-^±°-3 for a sample of 28 clusters without cooling flows. 
With cooling flows, a higher central value for fb and a shallower profile are reported, as would be 
expected. The fractional slopes of both these fits are steeper than that of the mean T profile of 
all the models. This is likely due to the lack of cooling in the simulations, which would deposit a 
concentrated amount of gas near the cores of the clusters. This gas would represent a small fraction 
of the total mass of the halo gas, but would elevate the baryon fraction in the interior. This should 
not greatly affect the results at radii near the virial radius, since the cooled gas would constitute a 
small fraction of the total gas mass at that point. 

These results indicate that the observed baryon fractions in galaxy clusters lead to under- 
estimates of the universal baryon fraction by a factor of 5% to 20% and that this factor is dependent 
on the cosmology. The anti-biasing of the gas is greatest in the hierarchical run, consistent with the 
picture of Pearce et al. (1994) in which the bias is due to transfer of energy from the dark matter 
to the gas during merging. Since cooling flows are not created in these simulations, the true bias 
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factor is not well established. However, cooling flows will only alter the baryon fractions in the 
inner radii since the total mass cooled is only on the order of 10^^ Mq or less. Ultimately, this result 
exacerbates the baryon overdensity problem in galaxy clusters. The variation in the enrichment of 
baryons with distance from the cluster centres also makes measurements of /(, problematic. It is 
crucial that the cluster be observed to radii of up to two virial radii to reduce the error to less than 
10%. 

Observations of the baryon fraction in clusters will be hampered by this biased distribution of 
gas in clusters. Measurements must be made to more than twice the virial radius, sampling regions 
for which accurate densities of both mass and gas are difficult. This uneven distribution is not as 
significant if there is any degree of smoothing, implying that any attempt to formulate a correction 
factor will be dependent on the cosmology at the 10% level. Fortunately, this is at the present level 
of measurement error. 
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Smoothing method 


f smooth 


Renormahsed 




(/i-^Mpc) 


to CTg 


unsmoothed 





N.A. 


3Mpc top-hat 


3 


Yes 


(7Mpc)-i k-cutoff 


7 


Yes 


7Mpc top-hat 


7 


No 


(14Mpc)-i k-cutoff 


14 


No 



Table 1: Properties of the simulations. Given are the effective smoothing radius and whether the 
data, after smoothing, were renormalised to erg = 0.935. 





-^clusters 


Mmin[lO'''MQ] 




unsmoothed 


100 


0.05 


1.52 


3Mpc top-hat 


48 


0.05 


2.85 


(7Mpc)-i k-cutoff 


48 


0.05 


1.75 


7Mpc top-hat 


8 


0.03 


1.95 


(14Mpc)-i k-cutoff 


6 


0.15 


2.27 



Tabic 2: Results of the cluster search. Given for each run is the number of clusters found, Ndustersi 
and the range of masses of the clusters, Mmin and M^ax- 
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Fig. 1. — The normalised baryon fraction as a function of mass. The baryon fraction is measured 
by T200 ^ a function of M2oo- 
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Fig. 2. — The normalised mean baryon fraction profiles, as measured by T(r/i?20o), for the models. 
At the top are the profiles for (T(r/i?20o)) at the specified radii. The mean cumulative enrichment 
profiles are given in the lower panel. 



